Dilepton low pr suppression as an evidence of the Color Glass Condensate 
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The dilepton production is investigated in proton-nucleus collisions in the forward region using 
the Color Glass Condensate approach. The transverse momentum distribution [pr), more precisely 
the low pt region where the saturation effects are expected to increase, is analyzed. The ratio 
between proton-nucleus and proton-proton differential cross section for RHIC and LHC energies is 
evaluated, showing the effects of saturation at small pr, and presenting a suppression of the Cronin 
type peak at moderate pt- These features indicate the dilepton as a most suitable probe to study 
the properties of the saturated regime and the Cronin effect. 

PACS numbers: 11.15.Kc, 24.85.-|-p 



I. INTRODUCTION 



At high energies, the Hnear evolution equations, based 
in the standard perturbative QCD, predict a high gluon 
density, requiring that the growth of the parton den- 
sity has to have a limit jlj - otherwise violation of the 
Froissart-Martin bound occurs - and expected to sat- 
urate at a scale Qs, forming a color glass condensate 
(CGC)ll H S 0. In this context, the search for sig- 
natures of a CGC description of the saturated regime 
is an outstanding aspect of investigation in heavy ion 
colliders. The first results of Relativistic Heavy Ion Col- 
lider (RHIC) , on charged hadron multiplicity in Au — Au 
collisions, were treated considering that the CGC for- 
mulation gives a natural qualitative explanation of the 
data'6'|. However, there are several issues to be clari- 
fied, before conclude that the dynamical of the partonic 
system should be described by a CGC already at RHIC 
energies. Particularly the current data are reasonably de- 
scribed by other models based on different assumptions 
"s"]. However, the charged multiplicity distribution in 
pseudo-rapidity for deuteron-gold collision is estimated 
within the CGC formalism at the deuteron fragmenta- 
tion region, and pointed out as a probable signature of 
the saturated regi me |9|1 . For a review of the CGC signa- 
tures see Refs. [ifl Ell- 
in order to investigate the high energy limit of the 
partonic interactions, the proton-nucleus scattering was 
proposed as an ideal experiment to give evidences of 
the saturation effects described by the CGC in the pro- 
ton fragmentation region 0,0,^3- Furthermore, the 
dilepton production was shown to be a sensitive probe of 
the perturbative shadowing and saturation dynamics in 
proton- proton, proton- nucleus and nucleus-nucleus scat- 
tering ill, lii 0, E 0, H 0i HI m in the forward 
kinematical region. It is an interesting observable since it 
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is a clean process in which there is no strong interaction 
with the nuclear medium final state. 

In this work we investigate quantitative features of 
the dilepton production in the forward region of proton- 
nucleus collisions in the context of the color glass con- 
densate. In particular, the transverse momentum (pt) 
distribution is studied focusing attention to the small px 
region, where the saturation effects are expected to be 
more significant. The main goal of this work is to show 
the effects of saturation at small px described by the 
CGC. Their presence at small pr should be considered 
as a possible signature of the saturation effects when con- 
trasted with proton-proton results. This comparison is 
performed evaluating the ratio between proton-nucleus 
and proton-proton cross section. This ratio shows two 
diferent behaviors; it presents Cronin type peak (if a lo- 
cal Gaussian for the correlator function is used) and a 
large suppression (if a non-local Gaussian is used), being 
a most suitable probe of the status of the Cronin effect 
as a final or initial state effect. This work is organized as 
follows. In the next section one presents a brief discus- 
sion on Color Glass Condensate formalism. In Sec. 3 the 
dilepton production cross section within the CGC formal- 
ism is presented. The Sec. 4 is devoted to the study on 
the color field correlator, which is a fundamental factor 
in the CGC approach. The numerical results are given 
and discussed in the last section where our conclusions 
are also presented. 



II. THE COLOR GLASS CONDENSATE 

The Color Glass Condensate picture holds in a frame 
in which the hadron propagates at the speed of light and, 
by Lorcntz contraction, appears as an infinitesimally thin 
two-dimensional sheet located at the light cone. The for- 
malism supporting this picture is in terms of a classical 
effective theory valid at small x region (large gluon den- 
sity) , and was originally proposed to describe the gluons 
in large nuclei 

At small X and/or large A one expects the transi- 
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tion of the regime described by the standard pertur- 
bative QCD (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP), BaUtsky-Fadin-Kuraev-Lipatov (BFKL)) to 
a new regime where the processes hke recombination of 
partons should be important in the parton cascade 0. 
In this regime, the growth of the parton distribution is 
expected to saturate below a specific scale Qs, forming 
a Color Glass Condensate ^ Qk i^, This saturated 
field, meaning the dominant field or gluons, has a large 
occupation number and allows the use of semi-classical 
methods. These methods provide the description of the 
small X gluons as being radiated from fast moving color 
sources (parton with higher values of x), being described 
by a color source density pa, with internal dynamics 
frozen by Lorentz time dilatation, thus forming a color 
glass. The small x gluons saturate at a value of order 
xG{x,Q'^) ^ l/as >> 1 for < Q^, corresponding to 
a multi-particle Bosc condensate state. The color fields 
are driven by the classical Yang-Mills equation of motion 
with the sources given by the large x partons. The large 
X partons move nearly at the speed of light in the positive 
z direction. 

In the CGC approach the light cone variables are em- 
ployed, where, = {x^,x~,x±), with x^ = l/^/2{x'^ ± 
x^) and x± = {x^, x^), and x^pp_ = x^p^+x^p^—x±-p±. 
The variable x~^ is the time light cone, and p^ its variable 

conjugated identified with the energy as p^ = ''"\p+^^ ■ 
The large x partons (fast) have momentum emitting 
(or absorbing) soft gluons with momentum A:+ << p+, 
generating a color current only with the -I- component 
= 5(x~)pa- In this framework, the nucleus is situated 
at x~ ~ 0, with an uncertainty Ax" ^ 1/^^, and there is 
a separation between fast and soft partons, implying that 
the former have large lifetime while soft partons have a 
short lifetime. These features assure that the color source 
density pa should be considered time independent, since 
for the emitted soft gluons (small x gluons) the source 
is frozen in time. However, after a time interval of order 
l/sp {Sp being the energy of the source in the light-cone) 
the configuration of the source is different. In order to 
have a gauge-invariant formulation, the source Pa must 
be treated as a stochastic variable with zero expectation 
value. For these reason, an average over all configurations 
is required and it is performed through a weight function 
Wa+ [p] , which depends upon the dynamics of the fast 
modes, and upon the intermediate scale A+, which de- 
fines fast (p"*" > A"*") and soft (p+ < A+) modes. The 
classical fields obey the Yang-Mills equation of motion. 



D,Frix) = 6^'+Paix-,x^), 



(1) 



and a physical observable is obtained by averaging the 
solution to this equation over all configurations of pa, 
with the gauge- invariant weight function Wa+ [p] ■ 

The effective theory is valid only at soft momenta of 
order A"*". Indeed, going to a much softer scale, there 
are large radiative corrections which invalidate the clas- 
sical approximation. The modifications to the effective 
classical theory is governed by a functional, nonlinear. 



evolution equation, originally derived by Jalilian-Marian, 
Kovner, Leonidov and Weigert ( JIMWLK) H |3| for the 
statistical weight function W\+ [p] associated with the 
random variable Paix). 

The solution for such functional evolution equation is 
not well determined yet and in order to make predictions 
or comparison with data, some phenomenological treat- 
ment should be given to the weight function. In this 
work an approximation to the weight function which is 
reasonable when we have large nuclei is used and con- 
sists in taking a Gaussian form [milll. As a con- 
sequence, most calculations in the CGC should be done 
quasi-analytically. In Ref. it is shown that a Gaus- 
sian weight function can accommodate both the BFKL 
evolution of the gluon distribution at high transverse 
momenta, and the gluon saturation phenomenon at low 
transverse momenta. A non-local Gaussian distribution 
of color sources has been predicted in Ref. l26j as a mean- 
field asymptotic solution for the JIMWLK equation and 
provides some modifications concerning phenomenologi- 
cal properties of the observables . The local Gaussian 
weight function assures that the color sources are corre- 
lated locally, on the other side, the non-local Gaussian 
allows correlations over large distances. In the following 
sections the phenomenological consequences of the choice 
of a local or non-local Gaussian type for the weight func- 
tion are investigated concerning the transverse momen- 
tum of the dileptons. 



III. DILEPTON PRODUCTION IN THE CGC 
APPROACH 

At high energies, the dilepton production in hadronic 
collisions looks like a bremsstrahlung of a virtual photon 
with momentum p decaying into a lepton pair, which can 
occur before and after the interaction of the quark (mo- 
mentum k) with the dense saturated gluonic field (mo- 
mentum q) of the target, in our case the nucleus A. We 
consider only the diagrams where the photon emission 
occurs before and after before the interaction with the 
nucleus, since the emission considering before and after 
the interaction vanishes '28'|. The dilepton production 
can be summarized in Figure ^ JJ„ JJ^ 21^ i 



9(k) 



9(q) 



Nucleus 

FIG. 1: Dilepton production in the CGC. 
Considering the Fig. the differential cross section 
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for the dilepton production in the CGC approach, for a colhnear quark {kx = 0), is written as |2ll |. 
I 



da. 



qA^ql + rX 



incl 



dzcPpTdlogM"^ 



{2n) 



[pI + M2(l - z)][{pT - zItY + ^'(1 - 



z(l - z)]Vf 



[pI + M^{\-z)] [{pT- zlTf+NP{\- z) 



(2) 



where fq represents the fraction of the electron charge 
carried by the quark q. The squared quark charge is 
— f^e^ and the charge from was incorporated 
in the in the expression. Ra is the nuclear radius, 
z = p~ /k~ is the energy fraction of the proton carried 
by the virtual photon, pT and are the transverse 
momentum and the squared invariant mass of the lepton 
pair, respectively; It = qr + Pt is the total transverse 
momentum transfer between the nucleus and the quark. 
The function C{It) is the field correlator function and 
defined by HI, 



C{It) = Jd 



(C/(0)C/t(x^))p, 



(3) 



with the averaged term representing the average over all 
configurations of the color fields in the nucleus, U{x± ) is 
a matrix in the SU{N) fundamental representation which 
represents the interactions of the quark with the classi- 
cal color field. The correlator considers the two diagrams, 
being the interaction at two transverse locations, and all 
the information about the nature of the medium crossed 
by the quark is contained in the function C{It)- In par- 
ticular, it determines the dependence on the saturation 
scale Qs (and on energy). 

In order to obtain a hadronic cross section, the va- 
lidity of the colhnear factorization in the fragmentation 
region is assumed and the expression in Eq. ^ is convo- 
luted with the partonic distribution function in the pro- 
ton (deuteron or nucleus) , as was done in ^3, and 
the cross section reads as, 



dp"^ dM dx F 



l^n2 ^ 



X2 



dli 



(27r)3 



ItW{pt,It,xi)C{It,X2,A), (4) 



where a;_F is the longitudinal momentum fraction given by 
xp — xi — X2, and xi and X2 are the momentum fraction 
carried by the quark from the proton and by the gluonic 
field from the nucleus, respectively. The expression Q 
is valid in the forward region, which means positive xp, 
or positive rapidities y, and the variables xi and X2 are 



defined by. 




or 



(5) 



(6) 



where M|. 



AP -\- p^ is the squared dilepton 
transverse mass and s is the squared center of 
mass energy. Here, using the structure function 
F2{x,Q^) = Y^z^^lMl^i^'Q"^) + qi{x,Q^)], the weight 
fimction W{pt,Iti Xi) can be written as. 



W{ptJt,xi)= dzzF2{xi/z,M^) 

J Xi 



(I-f (1 - Z)2)z2|2 



[p| + M2(l 

z(l - z)M^ 



zmpT-zlT)^ + A'P{l~z)] 



[p^^ + AP{l-z)] 



1 



[{pj.-zlj,)2+M2^l_,) 



(7) 



In our calculations the CTEQ6L parametrization [23 was 
used for the structure function, and the lepton pair mass 
gives the scale for the projectile quark distribution. The 
function W(j)t, It, xi) plays the role of a weight function, 
selecting the regions of dominance on It contributing to 
the cross section. 

In Eq. 0) the correlator function appears with an 
energy dependence (dependence on X2), which is not 
included in the original McLerran-Venugopalan model. 
One includes such dependence in the field correlator func- 
tion only in the saturation scale Qs,A — > Qs,Aix) to sim- 
ulate a low X evolution, as was done in the Ref. |23|, in 
order to investigate the effects of the x evolution in the 
dilepton pt spectra. The x dependence is parametrized 
in the form proposed by Golec-Biernat and Wiisthoff 
(GBW) ^ {Ql = (xo/x)^), with the parameters taken 
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from the dipole cross section extracted from the fit pro- 
cedm'e by GBW and CGCfit [sJl parametrizations, 
which will be discussed later. 
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FIG. 2: Weight function for lepton pair mass M = 3 GeV 
and rapidity y = 2.2 versus It- 



In Fig. [3 we plot the weight function for an specified 
lepton pair mass Af = 3 GeV, in the forward region (pos- 
itive xp\ with a positive value of the rapidity y = 2.2, 
considering the center of mass energy y^s — 350 GeV 
(RHIC). In such figure the results for three representa- 
tive values of px are presented, where a peak at It ~ Pt 
and a suppression at It < Pt are in order. These charac- 
teristics assure that the weight function selects the values 
of It larger than pT- Moreover, one verifies that larger 
values of pT provide a reduction in the normalization of 
the weight function at large values of It, when compared 
with the normalization at j3t = GeV. This pt behavior 
of the weight function is essential in order to determine 
the spectrum. As will be verified in the Sec. the pT 
distribution is suppressed for large values. 

All high density effects on the nucleus are encoded in 
the field correlator function. It is well known that the 
saturation effects in the correlator function C{IttX,A) 
are present below the saturation scale, meaning the small 
It region (In the next section one observes this feature in 
the Fig. 121 . Such behavior determines that only at small 
Pt the effects of saturation in the function C{ItiX,A) 
should be measurable, once the weight function selects 
the values of It larger than pT- 

To make a quantitative prediction to the dilepton pro- 
duction, the correlator function C{IttX,A) has to be de- 
termined. It plays an important role in the Color Glass 
Condensate formalism and should be compared with the 
dipole cross section. It is related to dipole Fourier trans- 
form and in the following, the function C{It,X2t A) will 
be discussed, through the analyses of some phenomeno- 
logical models. 



IV. THE COLOR FIELD CORRELATOR 

C{It,x,A) 

The function C{It) is considered a fundamental quan- 
tity in the CGC formalism, since it contains all the in- 
formation on high density effects. It can be related to 
the Fourier transform (F.T.) of the dipole cross section 
in the following way [22, H^H^ , 



C(Zt) = 



cPx^e''''^'''^[(Jdip{xi_ ^ oo) - adip{x^)],{?>) 



where ctq is the normalization of the dipole cross section 
at the saturation region {x± — > oo). Considering the 
GBW model for the dipole cross section adip[xji^,x) — 
CTo[l — exp(— (5^{x]x^/4)] [iOl the correlator function can 
be written as 



C{It,x,A) 



47r 



GBW 



Ql(.x,A) 



(9) 



where a simple dependence on energy {x) and atomic 
number {A) was taken into the saturation scale. Namely, 
the nuclear saturation scale was considered as Q^s{x, A) = 
A^/^Q^g{x) with Q^g{x) being the proton saturation scale 
parametrized of the form proposed by GBW Ql — 
{xo/x)^ I^Ol) where the parameters xq = 3.10"'* and 
A = 0.288 were determined from the fit to the Hadron 
Electron Ring Accelerator (HERA) data. This ansatz to 
the nuclear dependence of the saturation scale was stud- 
ied in the Ref. [S^l concerning the eA data, and was 
shown to be a consistent approximation for large nuclei 
and moderate energies. 

However, the GBW F.T. (Eq. ©) does not recover 
the perturbative behavior at large pt, since it presents 
an exponential tail, as we show in the Fig. |31 

Considering the McLerran-Venugopalan (MV) model, 
the function C{It) has no energy dependence and should 
be computed by taking the MV dipole cross section 2^, 



,(-*/ 



^(l-Jo(prj 



(10) 



The Fourier transform can be numerically computed in 
the form HI, 



C{It) 



MV 



^ d 



i^Xie^'-^-e-^/f^^'-""^^""", (11) 



where the value of is fixed depending on the energy. 
However, no x evolution is presented in the MV model 
and the energy dependence is introduced in the correlator 
by fixing the value of the saturation scale. Following the 
equation Hll|l . we propose to introduce the dependence 
on the energy and nuclei into the saturation scale in the 
form, 



C^MV^.ai^TjX^A) 



d x±e 



X e 
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The nuclear saturation scale is parametrized in the form 
presented previously, where the x dependence in the sat- 
uration scale (x) is considered from the parameters ex- 
tracted from the fit to HERA data and should be taken 
from GB W saturation model js^l , or from a dipole cross 
section based on the CGC approach (sj. Here it should 
be interesting to point out that in the recent work of Ref. 
[2^ , the Cronin effect was studied in the MV model, and 
the same energy dependence for the saturation scale was 
considered. Such a simple inclusion of the quantum cor- 
rections results in a disagreement with the RHIC d — Au 
data at forward rapidities concerning the Cronin effect 
[35I Is^ . In that work [23| , a non-local Gaussian distri- 
bution for (C7t(0)C/(a:±)) was introduced and the shape 
of the curves agrees with the Broad Range Hadron Mag- 
netic Spectrometer (BRAHMS) data at large rapidities, 
however presents large suppression at central rapidity. 
This disagreement shows that the dynamics of the CGC 
is a subject which deserves more comprehensive studies. 
Here we point out that the dilepton transverse momen- 
tum analyzed with local and non-local Gaussian correla- 
tor is a good observable to investigated such dynamics. 

In the large It {It » Qs) limit, the correlator func- 
tion should recover the perturbative behavior (l/Z^), and 
considering the MV model the correlator function can 
be expanded and written in a simple analytic expression 

CMV„„d(^T,a;, A)|;^>>Q^ ^2^-^^ — ^ 

I 



which emphasizes the large saturation effects in the small 
It region, as in the Fig. ((SJ. 

In a recent work, Ref. [sj it has been analyzed the 
structure function F2{x,Q'^) for x < 10~^ and 0.045 < 
< 45 GeV^, within the dipole picture, taking an ex- 
pression to the dipole cross section which interpolates the 
BFKL solution at r << 1/Qs{x) and the saturated be- 
havior at r >> 1/Qs{x), where the scattering amplitude 
saturates at one. The parametrized dipole cross section 
can be written as crdipix,r) = 2TiR^M{rQs,x) with [sij. 



A/'(rQ„y)=AAo(^) ^ ' to rQ, <2(^4) 

U{rQs,Y) ^ 1 - e-<^y^^\brQ,) > 2, 



where Y = ln(l/.T). There arc three free parameters: the 
proton radius R, the value xq of x at which the satura- 
tion scale has to be equal to 1, and the parameter which 
controls the energy dependence of the saturation scale A. 
The parameters a and b are determined to assure that 
Af is continuous at rQs = 2 (at least at first derivative). 
From the fit to the HERA data for the inclusive structure 
function F2{x,Q^) the parameters depend on the quark 
mass niq. 

Following this dipole cross section, we construct a func- 
tion C{It,X2, A) taking the Eq. ©, that will be called 
CGCfit, and obtain the following expression. 



CilT,x,A)cGC = 27r 



2/Q. 



rdrJo (hr) ( 1 - Wo cxp <{ 2 In ( ^ 



7.+ 



ln(2/rQ,) 



kA In 1 /x 



rdrJo(?T?-)e-"'"'^^'"^=) 



(15) 



2/Q. 



The energy and nuclear dependences are introduced with 

Qlia:,A)^A^/^f)\ 

Considering the two models for the dipole cross section 
(GBW and CGCfit) there is a set of parameters which de- 
termine the saturation scale, where the ones used in this 
work are presented in the Tablen]( the set of parameters 
are identified as fitl, fit2 and fit3), where the value of the 
saturation scale was calculated at a; = 10~^ for gold. It is 



shown that the recent CGC fit parameters provide a small 
value for the saturation scale, and this behavior should 
affect the dilepton production, as we will see in the next 
section. The nuclear radius is taken from the Woods- 
Saxon parametrization of the form, Ra = 1.2A^/^ fm, 
while the proton radius is taken from the fits and pre- 
sented in the Tabled 



In the Figure|31the function C{It,X2,A) is shown for 
the models discussed here and one verifies the large satu- 



ration effects at small It when we compare the functions 
obtained from the GBW and CGCfit with the asymp- 
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Parameter 


GBW 
(fitl) 


CGC fit m, = 10 MeV 
(fit2) 


CGC fit = 140 MeV 
(fit3) 


Xo 


3 X 10"'' 


1.06 X 10"" 


0.267 X 10"* 


X 


0.288 


0.285 


0.253 


Ql (x = 10-^ A = 197) 


4.114 GeV^ 


3.069 GeV^ 


2.327 GeV^ 


-Rp (Proton radius) 


0.6055 fm 


0.566 fm 


0.641 fm 



TABLE I: Parameters of saturation scale from GBW and CGCfit. 



totic behavior of the correlator function. At this point 
it is interesting to emphasize some features of the func- 
tions C{It -.x-i,^ A) extracted from the Fourier transform 
of GBW dipole cross section and the Fourier transform 
of the CGCfit. Considering the GBW F.T. (dashed-hne) 
the function C{It-,Xi,A)gbw depends on e~'r and is 
suppressed at large It (this behavior is shown in the inner 
plot on Fig. 13)1 . It results in an unrealistic suppression of 
the observable cross section at large , as emphasized in 
the Refs. [S^jll^- Considering the results from CGC fit 
(solid-hne), the function C{],t^xi-,-A) presents negative 
values at moderate 1^ (as can be seen in the Figure 
and this behavior should be due to the continuity only at 
first derivative or by the approximations in the construc- 
tion of the dipole cross section model 'si'l . Having those 
aspects in mind, the function C(It tX2, A) based on the 
McLerran-Venugopalan model, including an energy and 
nuclear dependence in the saturation scale, is employed 
here. The parameters to the saturation scale are taken 
from the fitl (triangle- up- line) and fit2 (circle-line) and 
one verifies in the Fig. O that the value of the satura- 
tion scale provides a small difference in the correlator 
obtained from the MS mod model at small It- 



C(It,x,A) {F,T, fit CGC dip) 

C(It,x,A) (F.T. GBW) 

. C(l„x,A) (iVIV„„ Q, CGCfit m,=10 iMeV) 
. C(i„x,A) (iVIV„„, Q. GBW) 
> C(i„x,A) (iVlV„., asymptotic) 




FIG. 3: Correlator C(It ,x, A) as a function of It- 

The correlator function is suppressed at large It-, and 
the weight function suppress the values of It smaller than 
Pt, the behavior of the cross section coming from the bal- 



ance between these two quantities. In such balance, the 
small pt dileptons clearly are the dominant contribution 
for the cross section and provide a physical probe for 
the CGC and for the models to the color field correlator 
function. 

The field correlator function presented up to here in 
MV model is obtained considering a local Gaussian func- 
tion for the weight function Wa+ \p\ . The correlator func- 
tion is defined by 



d^x I 



{U{0)U\x^)), 



(16) 



and the local Gaussian enters in the calculation of the av- 
eraged term (C/(0)C/^(xj^))p. The consideration of a non- 
local Gaussian function modifies the correlator in such 
way that it is written as 24, 27] 



C{It,x,A) = / (Pxx_ 



(17) 



with 



x{x,x±,A) = 



- I ^{i-Mx^p)) 

icj p 
ln(l+ f^li^)'^' 



(18) 



where, 7 is the anomalous dimension (7 w 0.64 for 
BFKL) and c « 4.84 This non-local field correla- 

tor function is plotted in the Figure0]in contrast with the 
same correlator obtained with the local Gaussian weight 
functional. 

The physical effect of the nonlocal Gaussian weight 
function is that the gluon sources are no longer corre- 
lated locally, as in the local Gaussian, but correlate over 
larger distances. This implies in a more drastic reduc- 
tion of the gluon density, as can be seen at small It in 
the Fig. 01 where the solid line represents the correlator 
function with a local Gaussian weight function and the 
long-dashed line represents the non-local Gaussian weight 
function. The effect of the local or non-local Gaussian 
weight function in the pt dilepton spectra will be dis- 
cussed in the next section in the context of the defined 
ratio RpA- 

Having addressed all the fundamental aspects to de- 
velop the calculation of the dilepton transverse momen- 
tum in the CGC formalism, one presents in the next sec- 
tion the numerical predictions using such approach and 
the discussions. 
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It (GeV) 

FIG. 4: Field correlator function with local and non-local 
Gaussian weight function. 



Figs. OandlBJ where the hne-diamond represents the cal- 
culation with the MV model, fixing the saturation scale 
at a value ^ 3 2 GeV ^ and = 8 GeV^, respec- 
tively. The X evolution provides a small suppression of 
the large px dilepton, in both cases. 
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V. RESULTS AND DISCUSSIONS 

In what follows, the numerical results on the dilep- 
ton transverse momentum distribution in CGC are ad- 
dressed and discussed. We consider pA collisions at 
RHIC (yi = 350 GeV) and LHC energies (^s = 8.8 
TeV) in the proton fragmentation region (positive ra- 
pidities). The calculations are performed fixing values of 
rapidities (or xp) and lepton pair mass M. We use the 
function C{It,X2) based on the McLerran-Venugopalan 
model, however an x dependence through the satura- 
tion scale is introduced, taking the parameters from the 
HERA data fit procedure GBW (fitl) [13 and CGCfit 
(fit2) (211 . For sake of comparison, the same differential 
cross section using the original McLerran-Venugopalan 
model, fixing a value to the saturation scale is evaluated. 

In Fig. 13 we present the transverse momentum dis- 
tribution for RHIC energies (-y/s = 350 GeV) in pA col- 
lisions, for a lepton pair mass M = 3 GeV and as in 
the Ref. at rapidity y = 2.2. The proton struc- 

ture function is taken from the CTEQ6L parametriza- 
tion [23. The solid line is the calculation with the 
McLerran-Venugopalan model, with the x dependence on 
the saturation scale, taking the parameters from the fit2; 
the dashed-line is the same calculation with the satura- 
tion scale taken from the fitl and the dot-dashed line 
is the calculation with the asymptotic behavior of the 
MV correlator function. Considering the transverse mo- 
mentum distribution at fixed mass and rapidities, the 
effects of quantum evolution are not too relevant in the 
range of transverse momentum investigated here, once 
the parametrization of the saturation scale assures that 
it is almost fixed in the region treated in this case, 
changing only weakly with the transverse momentum 



FIG. 5: Dilepton production at RHIC energies (y^ — 350 
GeV) in pA collisions, considering rapidity y = 2.2 and lepton 
pair mass M = 3 GeV. 

In the Figure|31the large saturation effects presented at 
Pt ^ 2 GeV are verified if one compares the asymptotic 
behavior of the correlator function with the MVmod pre- 
diction. As was shown in the last section, the asymptotic 
behavior of the correlator function {It >> Qs) depends 
on the Q'^/l^, then an increase of the saturation scale pro- 
vides an increase in the differential cross section at large 
Pt, as can be seen in the Figure [S] As a most interesting 
feature, only at large pt the effects of the choice of sat- 
uration scale affect the cross section, and the difference 
between the predictions being a factor of 2 considering 
the smallest value of the saturation scale, which is taken 
from the fitS, in contrast with the GBW ones. 

In Figure the dilepton transverse momentum distri- 
bution at LHC energies is shown, taking the same value 
of rapidity y ~ 2.2 to assure the forward region and to 
make a comparison with RHIC energies. The same be- 
havior concerning the saturation effects is verified, al- 
though such effects start to be significant for pt ^ 4 
GeV. The estimative with the MV model was performed 
and the suppression at large pt when the energy depen- 
dence is introduced in the saturation scale can be seen. 
The Pt spectra is enlarged at large pT if the saturation 
scale is large, as was verified for RHIC energies. 

In order to avoid any ambiguity with normalization, 
the ratio between the proton-nucleus and proton-proton 
differential cross section for RHIC and LHC is defined. 



daijpA) 



RpA = 



-KR^dMdxpdpl 



(19) 



{X2 



). Such behavior can be seen in the 



nR'^dMdxFdp^ 

Some attention should be given to the uncertainty in the 
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distribution for the weight function Wa+ [p] . For RHIC 
energies the solid hne represents the calculation for ra- 
pidity y — 2.2 and the dashed line for rapidity y = 3.2. 
For LHC energies the long-dashed line represents the cal- 
culation for rapidity y ~ 2.2 and the dot-dashed line for 
rapidity y — 3.2. It is verified that at moderate px the 
calculations show a Cronin type peak for RHIC and for 
LHC (there is a suppression of the dilepton production 
at RHIC and LHC energies comparing with the proton- 
proton collisions at small pt, at intermediate px the ratio 
is larger than 1, and there is a suppression at large px 
up to get the value 1). Such peak increases and is shifted 
to larger pT at larger rapidities, due to the fact that the 
saturation scale grows with the rapidity and no evolution 
is taking into account. 



FIG. 6: Dilepton production at LHC energies — 8.8 TeV) 
in pA collisions, considering rapidity y — 2.2 and lepton pair 
mass M = 3 GeV. 

determination of the nuclear radius, then each cross sec- 
tion is divided by the nuclear or proton radius. The factor 
A^/"^ was used in the denominator to guarantee a ratio 
RpA about 1 at large pr. 

The expression to the ratio RpA in the dilepton pro- 
duction defined here should be written of the form, 

RpA {y,PT) = 

J dPlrWipT, It, xi)Ca{It,X2, A) 

J d^lTW{pTjT,Xi)Cp{lT,X2)' ^ ' 

where Ca is the nuclear correlator function and Cp is the 
proton correlator function. The ratio in the Eq. (|20|l is 
similar to the one obtained in the Ref. |27[ to investigate 
the Cronin effect (Eq. (113) in the Refn27l|). 

The Cronin effect was discovered in the late's 70s [s^l 
and is related to the enhancement of the hadron trans- 
verse momentum spectra at moderated px (2-5 CeV) in 
comparison with the proton-proton collisions (the ratio 
between pA and pp present a peak at moderate pt)- The 
effect should be interpreted as being originated by the 
multiple scatterings of the partons from the proton prop- 
agating through the nucleus, resulting in a broadening of 
the transverse momentum of the initial partons. This in- 
dicates the Cronin effect as an initial state effect. The 
Cronin effect was measured by the RHIC experiments, in 
Au — Au and d — Au collisions, however, the theoretical 
approaches cannot describe the effect in all the range of 
rapidity measured by the collaborations 'ssj. Although 
the Cronin effect concerns the hadron transverse momen- 
tum spectra, it is also expected in the dilepton transverse 
momentum spectra, since the effect of multiple scatter- 
ings is an initial state effect. Moreover, the ratio obtained 
in this work is similear to that is used to investigate the 
Cronin effect ^7]. 

In the Fig. ||7J) one presents the results for the ratio 
RpA to RHIC and LHC energies considering a correlator 
field function C{It, x, A) obtained from a local Gaussian 
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FIG. 7: Ratio between proton-nucleus and proton-proton at 
RHIC and LHC energies by the CGC approach at distincts 
rapidities with local Gaussian distribution for the weight func- 
tion W\p\. 

Concerning the Cronin effect, such peak is enlarged at 
large rapidities if the local Gaussian correlator function 
with the same energy dependence implemented here is 
used ll^l, in complete disagreement with the BRAHMS 
experiment at forward rapidities [35l l36j . In the same 
ReL jl^l the Cronin effect was studied using a non- 
local Gaussian distribution for the weight function and 
the Cronin peak suppression is reached. However, there 
is a suppression on the normalization at centrality re- 
gion, which is not consistent with the RHIC data in cen- 
tral rapidity jJSj j?_6i , emphasizing that the the non-local 
Gaussian weight function should be the right physics at 
forward rapidities. 

On the dilepton side, the behavior of the ratio RpA, 
shows the same features of the Cronin peak at forward 
rapidities when investigated with the local Gaussian (the 
peak is enlarged and shifted to high pT at large rapidi- 
ties). However, the ratio RpA was also investigated with 
the correlator function obtained from a non-local Gaus- 
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sian weight function Wa+ [p], and is presented in the Fig. 
IHl The suppression of the ratio RpA is verified show- 
ing exactly the same features presented in the Cronin ef- 
fect I23, iseing a possible clean observable to study this 
property. Although, the Cronin effect was considered as 
a final state one in the Ref. in our analysis, the 

dilepton production seems to clarify this aspect. It was 
obtained that the Cronin type peak (or the suppression of 
the Cronin peak) in the dilepton px distribution appears 
as an initial state effect. In the Fig. ||SJ) the solid line 
represents the calculation for rapidity y — 2.2 and the 
dashed line for rapidity y = 3.2 at RHIC energies. For 
LHC energies the long-dashed line represents the calcu- 
lation for rapidity y = 2.2 and the dot-dashed line for 
rapidity y = 3.2. 

1.8 

. g ; RHIC (y=2.2) ; 

RHIC(y=3.2) 

LHC (y=2.2) 

LHC(y=3.2) 
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FIG. 8: Ratio between proton-nucleus and proton-proton at 
RHIC and LHC energies by the CGC approach at distincts 
rapidities with non-local Gaussian distribution for the weight 
function W[p]. 

At RHIC energies, the effect of suppression appears 
if the nonlocal Gaussian in the correlator color field is 
used, sugesting the measurement of such suppression, al- 
though the detectors should not be able to measure such 
behavior 40]. On the other hand, at LHC energies, the 
suppression of the ratio RpA reach large values of pt and 
such suppression increases with the rapidity. It is inter- 
esting to address that at LHC the experimental facilities 
provide a detection of dileptons in the forward region 
with transverse momentum above 1.5 GeV, depending 
on the rate of the signal of the observable and on the sig- 
nal from physics and machine sources '4l'|. This feature 
assures that at LHC energies such suppression behavior 



should be detectable too. 

VI. CONCLUSIONS 

In this work the large saturation effects described by 
the Color Glass Condensate in the dilepton production at 
small Pt region and the dependence of the large pt spec- 
tra on the saturation scale value are verified. Although 
at RHIC, the transverse momentum distribution of the 
dilepton should not be measurable in the very small pt, 
at intermediate px, the comparison between pp and pA 
cross section, provides a tool to study the Cronin effect 
and the dynamics of the Color Glass Condensate. 

Particularly, the dilepton transverse momentum distri- 
bution presents the suppression of the Cronin type peak, 
as observed in the inclusive observables, if a non-local 
Gaussian is used. Such behavior is observed in the Fig. 
©. At the LHC energies, at forward rapidities, the ef- 
fect of suppression increases (Fig. ©). Such large sup- 
pression at high energies gives an indication that dilepton 
transverse momentum provides a clear probe of the Color 
Glass Condensate description of the high energy hadronic 
interactions in the forward rapidity region. Moreover, 
the dilepton pT spectra should be used to investigate the 
properties of the Cronin effects and indicates it as an ini- 
tial state effect. Our results confirm the studies of Refs. 
[TtL ITsL l20j concerning the saturation effects. In addi- 
tion to this analysis, in a recent work fl^, the high pT 
and low mass region in the dilepton production in per- 
turbative QCD with all-order resummation was pointed 
as a good probe of the gluon distribution, as was indi- 
cated in the Refs. considering the dipole approach. 
Also, the mass distributions of the dileptons investigated 
in the CGC should identify effects of saturation at small 
mass region The ensemble of these features shows 
that dilepton production is an observable that deserves 
to be measured, once it carries information about the 
high density nuclear system. 
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